Statistical Estimation of the Quality of Quantum Tomography Protocols 
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I. INTRODUCTION 

Quantum information technologies rely on the use of 
quantum states in novel data transmission and comput- 
ing protocols [H-Q. Control is achieved by statistical 
methods via quantum state reconstruction. At present 
quantum state and process tomography serves as a prin- 
cipal instrument for characterization of quantum states 
preparation and transformation quality [4434j. 

In the last few years, a significant effort was made to 
improve tomographic methods of quantum state mea- 
surements. In [35( a method of choosing the best set 
of measurements from a finite list of possible ones is pre- 
sented. The basis of consideration for such problems is 
the asymptotic theory of statistical estimations of the 
parameters of the density matrix. Difficulties in solving 
optimization problems are caused by the large number of 
parameters to be estimated. In general case, it requires 
extensive numerical calculations. The number of possi- 
ble measurements becomes extremely large if one wants 
to cover a multidimensional parameter space of the mea- 
surement apparatus. Thus, holding a finite number of 
samples in multidimensional parameter space, does not 
guarantee that the truly optimal set of measurements 
is chosen. In other words, even the method is universal 
and works for qubits as well as for higher dimensional sys- 
tems, it does not say which particular measurements have 
to be performed. It just helps to choose the best mea- 
surements out of large suboptimal set of possible mea- 
surements. Another significant problem is the choice of 
an adequate parametrization for the quantum states (we 
will return to this issue in Sec. IVII) . 

The work [36[ provides a theoretical approach to op- 
timal measurements of qutrit systems, where optimal 
means the use of a minimal number of bases for mea- 
surements. As expected, the use of mutually unbiased 
bases is optimal in that sense, although its experimental 
realization, not discussed in the paper, might be far from 
the optimal one. 

Another recent work (37| addresses the problem of den- 



sity matrix reconstruction for large multiqubit states, 
where a complete state reconstruction is impossible due 
to the exponentially growing number of undefined param- 
eters with the increase of qubit number in the system. 
The solution given is in finding only the permutation- 
ally invariant part of the density operator, which gives a 
good approximation of the true quantum state in many 
relevant cases. 

In the present paper, we introduce a complete method- 
ology for statistical reconstruction of quantum states 
which is based on analysis of completeness, adequacy 
and accuracy of quantum measurement protocols [38l - 
l4lj | . Completeness of the quantum tomography protocol 
provides the possibility of reconstruction of an arbitrary 
quantum states (both pure and mixed) and it is assessed 
by means of singular value decomposition of a special 
matrix built upon operators of measurement. Consid- 
ered SVD allows us to introduce the condition number 
K, which characterizes the quality of quantum measure- 
ment protocol. Adequacy employs redundancy of a mea- 
surement protocol compared to the minimum number of 
measurements required for reconstruction. It is assessed 
by consistency of redundant statistical data and mathe- 
matical model, based on quantum theory. Accuracy of 
statistical reconstruction of quantum states is based on 
a universal statistical distribution proposed in (40| . 

It is worth mentioning that in real experiments the ac- 
curacy of the state reconstruction depends on two types 
of uncertainties: statistical and instrumental ones. If the 
total number of measurement outcomes (sample size or 
statistics) is large enough, the instrumental uncertain- 
ties dominate over fundamental statistical fluctuations 
caused_by the probabilistic nature of quantum phenom- 
Practically, the required statistics, allowing us 
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to exclude statistical fluctuations, depends on the tomo- 
graphic protocol itself and the total accumulating time 
needed for taking data. From this point of view it would 
be useful to point out simple and universal algorithms for 
the estimation of the chosen protocol on the design stage 
before doing experiments, as well as the sample size for 
desirable quality of the state reconstruction. 
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The method considered in the present paper has the 
following features: (1) It is well suited in the case of 
multi-qubit state tomography. (2) It accepts the recon- 
struction for mixed states of arbitrary rank as well as 
for pure states. (3) It allows to compare various quan- 
tum measurements protocols with each other and, more- 
over, with respect to the fundamental fidelity level. (4) 
It also motivates the experimenter to manage available 
resources in the best way as well as to choose an opti- 
mal quantum measurement protocol. The paper is orga- 
nized as follows. Section [TT] introduces the mathemati- 
cal apparatus for general quantum measurement proto- 
col (for states with discrete variables). In Sec. IIIII we 
discuss a generalized statistical distribution for fidelity 
by introducing specific random value which can be called 
the loss of fidelity. Important operational criterion and 
corresponding quantifier (so called conditional number) 
of the protocols' quality based on completeness are sug- 
gested. Section IIVI relates to analysis of specific exam- 
ples of tomographic protocols of qubits which are based 
on the polyhedrons geometry. Section IVl collects results 
of numerical and physical experiments for both pure and 
mixed states of single and pair of qubits. Since it repre- 
sents a paradigmatic example, here we restrict ourselves 
to polarization degrees of freedom only; however, any 
other degrees of freedom can be considered. In Sec. IVII 
we discuss in details the features and advantages of the 
suggested approach in context of completeness, adequacy 
and accuracy. Then we conclude with Sec. I VIII 



II. QUANTUM MEASUREMENT PROTOCOL 

An arbitrary s-dimensional quantum state is com- 
pletely described by a state vector in a s-dimensional 
Hilbcrt space when it is a pure state, or by a density ma- 
trix p for a mixed one. To measure the quantum state 
one needs to perform a set of projective measurements 
on a set of identical states. 

A quantum measurement protocol can be defined by a 
so-called instrumental matrix X that has m rows and s 
columns [l7l - ll~9l |. where s is the Hilbert space dimension 
and m the number of projections in such space. For every 
row, that is, for every projection, there is a corresponding 
amplitude Mj, 

Mj=X jlCh j = l,2,...,m (1) 

where we assume a summation by the joint index I, the 
ci(l = 1, 2, . . . , s) being the components of the state vec- 
tor in the Hilbcrt space of dimension s. The square of 
the absolute value of the amplitude defines the intensity 
of a process, which is the number of events in 1 s 

*j = \Mj\ 2 . (2) 

The number of registered events kj is a random variable 
Poissonially distributed, tj being the time of exposition 



of the selected row of the protocol and Xjtj the average 
value, 

P(%) = ^-Uxp(-A^). (3) 

It is convenient to introduce special observables Aj (so 
called intensity operators), which are measured by the 
protocol during experiment, 

A, = tr(A jP ) (4) 

Here Aj = XjXj is the intensity operator for the quan- 
tum process Xj (the row of the instrumental matrix X). 
In this case the intensity operator for quantum process 
Aj is a projector, so we have 

A| = Aj (5) 

Formally, in the most general case, Aj is an arbitrary 
positively defined operator [42| ■ It can be presented as a 
mixture of projection operators described above. 

A, = 5>*j fc)t *« (6) 

k 

Here the index k sums different components of the mix- 
ture that have weights fa > 0. Such measurement can be 
conveniently presented as a reduction of the set of projec- 
tion measurements where only total statistics is available, 
while statistical data for individual components are not 
available. The general projection measurement is a par- 
ticular case of Eq. (J6]) where fa = 1, fa = fa = ■ ■ ■ = 0. 
If the sum of the intensities multiplied by the exposition 
time is proportional to a unit matrix, then we say that 
the protocol is brought to a decomposition of unity [42j , 

m 

I = J2^3=IoE, (7) 

where Iq is the constant which defines overall intensity 
and E is the identity (or unit) matrix. A protocol for 
which the condition ([7]) holds in the general case can be 
brought to the so-called non-orthogonal decomposition 
of the unity [42|]. In this case the protocol analysis is 
simplified. In mathematics such measurements are con- 
sidered as the most general extension of traditional von 
Neumann measurements, which arc based on the orthog- 
onal decomposition. Even if it is reasonable to require 
([7| due to the total probability preservation, it is worth 
mentioning that real experimental protocols often cannot 
be brought to decomposition of unity. Indeed, in real ex- 
periments when event registration scheme is used, the 
experimenter often adjusts his device to distinguish only 
one projection of the quantum state (simply loosing data 
which corresponds to the other states). Thereby real ex- 
periments (due to technical requirements) usually do not 
provide registration of the whole statistical ensemble and 
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thus they are not restricted by total probability preser- 
vation requirement. However, the method suggested in 
the present paper is applicable to these cases as well. 

The normalization condition for the protocol defines 
the total expected number of events n summarized by all 
rows: 



3=1 



(8) 



where tj is the acquisition time. Condition ([5} substi- 
tutes the traditional normalization condition for the den- 
sity matrices, tr(p) = 1 

In the following we consider the protocols of quan- 
tum measurements in terms of two important notions- 
completeness and adequacy. For this purpose, we intro- 
duce some opportune notation. First, when describing 
the whole sequence of quantum protocol measurements 
each quantum process intensity matrix Aj of dimension 
s x s is pulled into a single string of the length s 2 (to 
put the second string to the right from the first etc). 
Then, we assign a weight defined by an exposition time 
tj to each row in Bj and we construct a matrix m x s 2 
from these rows, calling it the measurement matrix of the 
quantum protocol. We assume that m> s 2 . 

In the case of projection measurements defined by rows 
Xj (j — 1, . . . , m) of instrumental matrix X, the rows 
Bj of the measurement matrix B could be calculated 
through the use of the tensor product of the row Xj and 
the complex conjugate row X* . 



B, 



tj ■ Xj 



Xj. 



(9) 



In the following exposition times are assumed to be equal 
to 1. With this matrix B, the protocol can be compactly 
written in the matrix form: 



Bp = T 



(10) 



with p being the density matrix, given in the form of a 
column (second column lies below the first, etc.). The 
vector T of length m records the total number of reg- 
istered outcomes. The algorithm for solving Eq. (fTU|) 
is based on the so called singular value decomposition 
(SVD) [43}. SVD serves as a base for solving inverse 
problem by means of pseudo-inverse or Moore-Penrose 
inverse (43l . |44| . In summary, the matrix B can be de- 
composed as: 



B = UVS\ 



(11) 



where U, (m x m), and V, (s x s ), are unitary matri- 
ces and S, (m x s 2 ), is a diagonal, non- negative matrix, 
whose diagonal elements are "singular values" . Then (|T0|) 
transforms to a simple diagonal form: 



Sf = Q 



(12) 



with a new variable /, unitary related to p via / = V' p, 
and a new column Q, unitary related to the vector T 



by the equation Q = WT. This system is easy to solve 
because S is a diagonal matrix. Its analysis allows classi- 
fying measurements from the viewpoint of adequacy and 
completeness [41 1. 

Let m > s 2 , that is, the number of measurements is 
greater than the number of elements in the density ma- 
trix. The rank of the model q denotes the number of non- 
zero singular values of the matrix B. By defining q we 
formulate two important conditions of any tomography 



protocol, namely its completeness and adequacy [41| . It 
is obvious that q < s 2 . The last m — q rows in the matrix 
S are equal to zero. Then it follows that for the system 
to be adequate, it is necessary that the last m — q values 
in the characteristic column Q are also equal to zero. We 
will refer to this condition as a measurements adequacy 
condition. If it does not hold the model is inadequate, 
that is, statistical data do not correspond to any quan- 
tum mechanical density matrix. It may mean, for ex- 
ample, that either the experiment is realized incorrectly 
or measurement matrix is wrong. Adequacy means that 
the statistical data directly correspond to the physical 
density matrix (which has to be normalized, Hcrmitian 
and positive). However it is worth noticing that, gener- 
ally, for mixed state it can be tested only if the protocol 
consists of redundant measurements (i.e. if to > s 2 ). 

Suppose that the model is adequate. The protocol 
is supposed to be informational complete if the number 
of tomographically complementary projection measure- 
ments is equal to the number of parameters to be es- 
timated; mathematically completeness means q = s 2 . If 
all singular values are knowingly nonzero, that is, q = s 2 , 
then unconditional completeness holds and a solution ex- 
ists and is unique. Measurement protocol completely re- 
turns any quantum state (pure and mixed) that could be 
defined in the considered Hilbert space. 

In this case we can determine the factor column divid- 
ing the elements of characteristic column by the corre- 
sponding singular values 



fj = Qj/Sj j = 1,2,. 



(13) 



As a result we obtain the desired density matrix by a 
unitary transformation: 



P = Vf. 



(14) 



Due to the unitarity of the matrix V , the factor column 
/ determines the degree of purity of the quantum state: 



3=1 



(15) 



Finally, suppose q < s 2 , that is, some singular values are 
equal to zero. In this case for nonzero values we have: 

fj = Qj/Sj, J 1.2 q. (16) 

Let us call these factors fj, j = 1,2, ... ,q defined factors. 
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At the same time, for nonzero values, we have equa- 
tions corresponding to uncertainty "zero divided by zero" 



0/.- = 0, J=q + l,...,S 2 



(17) 



We shall call these factors /, j = q+ 1, . . . , s 2 as unde- 
fined factors. As solutions of the last equations arbitrary 
complex numbers could be used. The considered situa- 
tion corresponds to the incompleteness of measurements 
and this system of equations has an infinite number of 
solutions. However, not all of them correspond to real 
physical density matrices. The physical solutions only 
correspond to that having a Hermitian nonncgativc def- 
inite density matrix. Formally all these solutions could 
be obtained by scanning all possible values of undefined 
factors. It is evident that such a procedure could be car- 
ried out only when the dimension of undefined factors' 
space is relatively small. 

Let us call regularized (normal) a solution correspond- 
ing to the special choice of all undefined factors: fj = 
for j = q + 1, . . . , s 2 . Due to the unitarity of relation be- 
tween the density matrix and the factor column, a regu- 
larized solution corresponds to the minimum purity level 
of the reconstructed state. In this case we have: 



tr(p 2 )>X> 



(18) 



One could see that in the case of incomplete protocols any 
additional measurements could cither cause an increase 
of the purity or leave it unchanged. If the regularized 
solution already describes a pure state, then new mea- 
surements will not influence the reconstructed state. In 
other words in this case we can obtain complete infor- 
mation about the considered state despite of incomplete 
measurement protocol. We will call such protocols con- 
ditionally complete, that is, there is completeness under 
condition that only specially chosen states are considered; 
for instance, this is the case of experiments like "which 
way" [H. 

Then wc assume that there is unconditional complete- 
ness (q = s 2 ). Equations (fT5)l and (jT4"f could be used to 
approximate the reconstruction of the density matrix if 
we substitute experimental events (frequencies) into the 
right side of Eq. (|TU)) . However due to statistical fluctua- 
tions of experimental data the reconstructed matrix will 
not always be positive definite (components with small 
weights could be reconstructed as negative numbers in 
this case). Despite this disadvantage, the method pro- 
vides a good zero approximation for widely used maxi- 
mum likelihood method (ML). In this case, the compo- 
nents with negative weights simply assumed to be zero, 
then the density matrix is multiplied by a factor that en- 
sures the correct normalization. Note that ML method it- 
self is free from the considered disadvantage because pos- 
itive definiteness lies in the nature of the method. At the 
same time the zero approximation, obtained from pscudo 
inversion method, significantly accelerates the search of 
ML solution. 



III. UNIVERSAL STATISTICAL 
DISTRIBUTION FOR FIDELITY LOSSES AND 
THE MAXIMUM POSSIBLE FIDELITY OF 
QUANTUM STATES RECONSTRUCTION 

The accuracy of quantum tomography can be defined 
by a parameter called fidelity [l], |4g| 



F 



(19) 



where p^ is theoretical density matrix and p is recon- 
structed density matrix. The fidelity shows how close the 
reconstructed state is to the ideal theoretical state: the 
reconstruction is precise if the fidelity is equal to one. 

This equation looks quite complex, but it becomes sim- 
ple if we apply the Uhlmann theorem [46| . According to 
this theorem, the Fidelity is simply the maximum possi- 
ble squared absolute value of the inner product: 



F=\{c \c)Y 



(20) 



where cq and c are theoretical and reconstructed purified 
state vectors. 

We explicitly use the Uhlmann theorem in our al- 
gorithm of statistical reconstruction of quantum states, 
based on maximum likelihood method. This fact is very 
important even if the state is not pure, we have to purify 
it by moving into a space of higher dimension (40| . 

It is well known that purified state vectors are defined 
ambiguously. However, this ambiguity does not preclude 
from reconstructing a quantum state, which is a very 
useful feature of the suggested algorithm. It is devised 
in the way that different purified state vectors produce 
the same density matrix and therefore the same fidelity 
during the reconstruction. This is a key principle for 
the proposed procedure and thus reconstruction can be 
obtained by means of purification. Purification greatly 
facilitates the search of a solution, especially when we 
need to estimate a great number of parameters (hundreds 
or even thousands). 

The fidelity level ([20]) has a simple probabilistic inter- 
pretation. If we choose a known reconstructed vector 
and its orthogonal complement as a measurement basis 
for an unknown state cq then F returns the probability 
that this unknown state coincides with the reconstructed 
one. 

It is equally important that due to the usage of purifi- 
cation procedure we succeed in formulating a generalized 
statistical distribution for fidelity [4(J. Purification pro- 
cedure will be described below in Sec. I VII [formula ([55)1 
provides with a transition from the density matrix to a 
purified state vector] . The value 1 — F can be called the 
loss of fidelity. It is a random value and its asymptotical 
distribution can be presented in the following form: 



1 - F 



(21) 
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where dj > are non-negative coefficients, £j ~ iV(0, 1) 
j = l,...,j m ax ar e independent normally distributed 
random values with zero mean and variance equal to one, 
Jmax = (2s — r)r — 1 is the number of degrees of freedom 
of a quantum state and corresponding distribution, s is 
the Hilbert space dimension, and r is the rank of mixed 
state, which is the number of non-zero eigenvalues of the 
density matrix. In particular j max = 2,s — 2 for pure states 



and j r 



s — 1 for mixed states of full rank [r = s) 



This distribution is a natural generalization of the \ 2 
distribution. Ordinary % 2 distribution corresponds to the 
particular case when d\ = d% = ■ ■ ■ = dj max = 1 (all com- 
ponents of vector d are equal to one) . In the asymptotic 
limit considered by us, the parameters dj are inversely 
proportional to the sample size n, that is, dj ~ — . This 
dependence allows for an easy recalculation when we use 
different sample sizes. 

The method of calculating the parameter vector d is 
based on Fischer's information matrix. In this case, the 
eigenvectors of information matrix define directions of 
principal fluctuations of the purified state vector whereas 
the respective variances of the principal fluctuations are 
inversely proportional to the eigenvalues of the informa- 
tion matrix. The method is described in (40j in details. 
From Eq. ([2"T]l one gets the average fidelity loss 



<l-F) = Y^dj. 
j'=i 



(22) 



It is also easy to show that the variance for the fidelity 
loss is 



Let us consider a special case of the quantum state, 
which is defined by a uniform density matrix. Moreover, 
this matrix is proportional to the identity matrix. Such 
state can be represented as a "white noise" for which 
all weights of principal components are equal. Let us as- 
sume that the protocol can be brought to projection mea- 
surements, which form non-orthogonal decomposition of 
unity in accordance with §5§ and (0). In this case there 
is a simple relation between the vector d, of size s 2 — 1 
determining the distribution of losses and the vector of 
the singular values for the measurement matrix B (the 
element of the greatest value should be eliminated from 
this vector of size s 2 ). Denoting the reduced vector of 
the size s 2 — 1 as b, one gets the relation between these 
vectors: 



11 ' W i = 1 ' 2 " 



,s 2 -l, 



(28) 



where the constant C is given by 



C 



3 

i( s -iy 



(29) 



For multiqubit protocols considered in the next section, 
which are based on polyhedrons, we have: 



dj = 



Asnb 2 ' 



(30) 



where m is the polyhedron's faces number and I is a num- 
ber of qubits in the register. It is obvious from (|30l) that 



- 2 = 2E4 

3=1 



(23) 



Moments of higher order for this distribution can be cal- 
culated analytically. For example the momentum of third 
order is called skewness and describes, for a random vari- 
able x, the asymmetry 



Ih 



M{[x- M(x)} 3 } 



(24) 



M denoting the mathematical expectation. The fourth- 
order moment is called excess kurtosis, 



M{[x-M(x)n 
Pi = 1 4 



(25) 



The equations for f3\ , /?2 are: 
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j'=i 



48 £ 4 



(26) 



(27) 



(31) 



In addition, it can be shown that for such protocols the 
following equation holds: 



(32) 



In this case the condition number K of the matrix B is 

i 



K = cond(B) = ( y/S 



(33) 



Recall that the condition number of a matrix is the ra- 
tio of the maximum singular value to the minimum one. 
Note also that in the definition all singular values 
are taken into account while in the (|32|) the one is ignored 
due to normalization. 

Let us introduce a value of fidelity loss, which is inde- 
pendent on the sample size. 



L = n(l-F)=nJ2d r 

3=1 



(34) 



This quantity serves as the main figure of merit of the 
precision in the examples analyzed below: the lower the 
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value of the loss function (|34[) . the higher is the preci- 
sion of the protocol. As we can see from the Eq. (f3"4"|). 
this value is determined by vector d that defines the gen- 
eral fidelity distribution. Therefore, the general fidelity 
distribution serves as a tool for completely solving the 
problem of precision for quantum tomography. 

As an important example, let us consider the proto- 
col defined by projection measurements, which form non- 
orthogonal decomposition of the unity. It can be shown 
that the following condition holds in this case 



— V — = s - 1 

4n ^ dj 

j=i j 



(35) 



Here v is the number of parameters to define a state: 



(2s - r)r - 1 



(36) 



Optimization of the problem leads to search of the mini- 
mum level of fidelity losses (|2"Tj) such that ([33]) holds. Ob- 
viously the problem is solved when d\ = d% = • • • = d v , 
when mean losses approach their minimum level: 



4n(s- 1) 



(37) 



Note that the requirement d\ — d 2 = ■ ■ ■ = d v not only 
defines a minimum level of mean losses ([22)1 . but also 
a minimum of other moments of losses [variance (|23|) . 
skewness (|26jl and excess kurtosis (|27j) ]. 

It appears from (|37j) that the minimum possible loss is 
given by the following equation: 



- opt 



(38) 



mm 4(s-l) 
For pure states v = 2s — 2, the possible loss is given by: 

1. (39) 
= s' 2 — 1 the loss is given 



opt 
min 



For mixed states of full rank v 
by the following equation: 

L opt = (S+1) 2 (S~1) 



(40) 



Any protocol for any quantum state cannot have losses 
lower that those defined by this equation, if the protocol 
can be brought to projection measurements, which form 
non-orthogonal decomposition of the unity. Note that if 
the protocol cannot be brought to decomposition of unity 
then the losses can be lower than defined by this equa- 
tion. However, this is true only for certain states and an 
improvement in reconstruction precision for some states 
is completely compensated by a significant deterioration 
in reconstruction precision for other states. 



IV. PROTOCOLS BASED ON THE 
POLYHEDRONS GEOMETRY 



this section are formed by projective quantum measure- 
ments on states that are tensor products of singlc-qubit 
states. For example, if a singlc-qubit measurement pro- 
tocol is formed by a projection onto a polyhedron with 
m faces inscribed in the Bloch sphere, then it has m rows 
(see below). Therefore, the corresponding /-qubit proto- 
col possess m l rows. 

Intuitively, projection of the state under consideration 
onto symmetric solids inscribed in the Bloch sphere leads 
to better accuracy of the reconstruction. Shown below 
is a study based on the calculation of the loss function 
(1341), w hich significantly extends the results obtained in 
[23T . [26j | . In the term of single-qubit protocols it is worth 
highlighting regular polyhedrons and polyhedrons with 
lesser but still rather high level of symmetry. 

Regular polyhedrons or Platonic solids are used for the 
most symmetrical and uniform distribution of quantum 
states on the Bloch sphere. Projections are defined by 
directions from the center of Bloch sphere to centers of 
polyhedron faces. Therefore, the number of polyhedron's 
faces defines the number of protocol's rows and is equal 
to 4 for tetrahedron, 6 for cube, 8 for octahedron, 12 for 
dodecahedron and 20 for icosahedron. 

These five bodies form the complete set of regular poly- 
hedrons. The search of quantum measurements proto- 
cols with high symmetry on Bloch sphere and number of 
rows greater than twenty requires to consider non-regular 
polyhedrons, which have high symmetry. As examples 
of such polyhedrons we have chosen fullerene (truncated 
icosahedron) that defines quantum measurement proto- 
col with 32 rows (equal to the number of fullerene's faces) 
and also a dualQto fullerene polyhedron (pentakis dodec- 
ahedron) which defines quantum measurement protocol 
with 60 rows (that is the number of its faces and also the 
number of vertices of fullerene) . 

It is noteworthy that all protocols considered here can 
be brought to decomposition of the unity [36] . 

A comparison of the maximal possible fidelity with the 
fidelity of the protocols considered here shows that as the 
number of polyhedrons' faces increases fidelity rapidly 
converges to the theoretical limit (in addition, rapidly 
increases a uniformity of fidelity distribution on the Bloch 
sphere). We should mention that the accuracy of the 
suggested protocols is much higher in comparison with 
protocols exploiting not so highly symmetrical states. 

As an example we calculated numerically a set of pic- 
tures that demonstrates the Bloch sphere scanning by 
means of various measurement protocols for the single 
qubit pure states. The corresponding colors indicate the 
value of the Fidelity loss function. 

Figure [1] defines the value of the loss function for the 
protocol based on tetrahedron. The minimal losses are 
equal to 1 as well as for other figures. The maximum 
losses for tetrahedron are equal to 3/2. 



In this section we present analysis of specific examples 
of tomographic protocols which are based on the polyhe- 
drons geometry. The multiqubit protocols described in 



1 A pair of polyhedra are called dual, if the vertices of one corre- 
spond to the faces of the other. 
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Tetrahedron - 4 rows, Lmin=1, Lmax=3/2 




FIG. 1. (Color online) Tetrahedron and distributions fidelity 
loss over the Poincare-Bloch sphere for the protocol based 
on tetrahedron geometry. Color bar shows level of average 
fidelity loss L. 

Figure [5] presents a cube and an octahedron. These 
polyhedrons are dual to each other. The maximum losses 
are equal to 9/8 in both cases. 

Figures [3] and 2] present a dodecahedron and an icosa- 
hedron, as well as fullerene and a polyhedron that is dual 
to the latter: one can clearly see that when the number 
of projections grows the maximum losses converge to the 
minimum possible losses. In the limit of infinite number 
of points on the Bloch sphere we get an optimal pro- 
tocol for which the precision of reconstruction does not 
depend on the reconstructed state at all. The price for 
that would be an infinite growth of the performed mea- 
surements, so in a real experiment one should sacrifice 
the fidelity losses (accuracy) to use a realistic (limited) 
number of measurements. However, this point is quite 
common for any quantum tomography protocols exploit- 
ing redundant (with respect to dimension of the recon- 
structed state) number of measurements. 

For the sake of completeness, fullerene and its dual 
protocols are presented on the last picture, while Table U 



Cube(6rows) and octahedron(8rows) Lmin=1, Lmax=9/8 




FIG. 2. (Color online) Shapes of solids and distributions of 
average fidelity loss over the Poincare-Bloch sphere for the 
protocol based on cube and octahedron geometry. Color bar 
shows level of average fidelity loss L. 



Dodecahedron - 12 rows Lmin=1, Lmax=36/35 




Icosahedron - 20 rows Lmin=1, Lmax=45/44 




FIG. 3. (Color online) Shapes of solids and distributions of 
average fidelity loss over the Poincare-Bloch sphere for the 
protocol based on dodecahedron and icosahedron. Color bar 
shows level of average fidelity loss L. 



presents the results of numerical experiments for pure 
quantum states with the number of qubits from 1 to 3. 
It is worth noting that the algorithm of numerical opti- 
mization does not garantee the finding of the global op- 
tima. Using the numerical procedure we have only found 
hypothetical maximum loss values. 

The precision of each protocol can be characterized by 
the following bounds L m in < L < £ max . This inequal- 
ity defines a rather narrow range, where the precision of 
quantum state reconstruction is localized definitely. 



Fullerene - 32 rows Lmin=1 , Lmax=234/233 




FIG. 4. (Color online) Shapes of solids and distributions of 
average fidelity loss over the Poincare-Bloch sphere for the 
protocol based on fullerene and its dual polyhedron. Color 
bar shows level of average fidelity loss L. 
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Numerical calculations demonstrate that the minimum 
possible losses L m i n are defined by the theoretically de- 



rived optimal limit L n 



L 



Op t 



1. 



By looking to the upper limit L ma x, shown in the Ta- 
ble HI as the result of numerical experiments, one can 
evince that for single-qubit protocols as the number of 
projections grows the maximum losses converge to the 
minimum possible losses. In the limit of an infinite num- 
ber of points on the Bloch sphere one gets an optimal pro- 
tocol for which the precision of the reconstruction does 
not depend on the reconstructed state. 

However, this is not true for multiqubit protocols. For 
two-qubit protocols, the maximum possible losses ap- 
proach the level L max m 3.38 while the minimum pos- 
sible level is equal to L^ n = 3. Similarly for three-qubit 
states the values are L max m 7.7 and L^ n = 7. These 
results are due to the fact that these protocols are based 
on projections only onto non-entangled states. 

Here we should stress again that from the theoretical 
point of view in the multiqubit case the discussed proto- 
cols are not the best possible ones because they do not in- 
volve projections onto entangled states. In that case the 
precision will be somewhat smaller than the minimum 
possible limit. (L max ss 3.38 compared to = 3 for 

two-qubit states and L max « 7.7 compared to L°^ n = 7 
for three-qubit states). 

It is also worth noting that, though the protocols based 
on polyhedrons with small number of faces (tetrahedron, 
cube) are somewhat less precise, they are much easier in 
practical implementation. 

When considering tomography of mixed states it is 
worth noting that there is no finite upper limit for preci- 
sion losses (losses can be infinitely large L — > oo). Such 
large losses arc inherent to mixed states, which are close 
to pure ones. In fact the number of real parameters that 
define a mixed state of full rank in Hilbert space of di- 
mension s is equal to s — 1, which is significantly greater 
for large s than for a pure state that takes only 2s — 2 real 
parameters. In case of a mixed state that has one pre- 
dominant component, the smaller weight components al- 
most do not affect the statistical data and do not increase 
the amount of Fisher information for reconstructing the 
greatly larger number of parameters. Theoretical analy- 
sis, numerical and real physical experiments completely 
prove these statements. 

The lower limit for precision losses can be applied to 
mixed states as well. In this case optimal minimum losses 
are realized for "white noise states" (uniform density ma- 
trix), when all components have equal weights. Then, 
there is a simple relation between the vector d of dimen- 
sion s 2 — 1, that defines distribution of precision loss, and 
the vector of singular values of the measurement matrix 
B. Corresponding estimate for L m j n for the protocols 
considered here is given by the following equation: 



10' 



1 



This value depends on the number of qubits, but does 
not depend on the type of polyhedron. It defines the 
minimum possible losses for the considered protocols that 
do not use projections on entangled states. 

Recall that in the general case for any protocols includ- 
ing those ones that involve projections onto entangled 
states minimum (optimal) losses are described by: 



- opt 



(2 ( + l) 2 (2 ( -1) 



4(« - 1) 



(42) 



(41) 



By comparing the Eqs. (|4"Tj) and (|4"2")l one can see 
that these protocols provide minimum possible (optimal) 
losses for reconstruction of mixed states of full rank only 
for single-qubit states. Therefore, for multi-qubit cases 
protocols that provide minimum possible losses during 
quantum states reconstruction should necessarily include 
projections onto entangled states. 

Let us consider few examples demonstrating the fea- 
tures of the developed approach. Figure [5] presents the 
results of numerical experiments testing the universal 
statistical distribution for fidelity; 200 experiments were 
conducted with sample size 1 million each. 

The measurement protocol is based on tetrahedron. 
We considered a four-qubit state that represents a mix- 
ture of GHZ state and uniform density matrix (white 



p = /^ + (l-/)|GHZ)(GHZ| ! (43) 

where E is the unit matrix of size 16 x 16, |GHZ) 
is the state of Greenberger-Horne-Zeilinger: |GHZ) = 
-i=(|0000) + 1 1111)), and / is the weight of the uniform 
density matrix (white noise). In our case / = 0.5 (50%). 
It is a multiparametric distribution, being the size of the 
vector of parameters 255. These data (Fig. [5]) demon- 
strate a good agreement between results of a numerical 
experiment and the theory with high critical significance 
level (0.65) for x 2 criterion. 

Then we consider the dependence of reconstruction 
precision on the weight of the "white noise" component. 
The value of fidelity can belong to a wide interval, thus 
it is convenient to use a new variable z = — log(l — F). 
Here and below log denotes the common logarithm. The 
new variable z defines the number of nines in numerical 
representation of fidelity, for example, z = 3 means that 
F = 0.999. 

As an another example, Fig. [5] shows distributions cal- 
culated by using the new variable for a three-qubit state, 
which is again a mixture of GHZ state and "white noise" 
(uniform density matrix). Calculations were performed 
with a protocol based on dodecahedron. The sample size 
n is equal to one million as well. It is evident that the 
higher the white noise weight the higher the precision of 
reconstruction. It is not difficult to explain this fact in 
the context of previous discussion. Namely the "white 
noise" is the best one for mixed state reconstruction. 
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TABLE I. Results of numerical experiments that define maximum precision losses L D 
geometry. 
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FIG. 5. (Color online) The universal statistical fidelity loss 
distribution for 200 numerical experiments (histogram) com- 
pared with the theoretical expectation (solid curve). Sample 
size of each experiment is 1 million. 




8 qubits (510 parameters ) 

7 qubits (254 parameters ) 

qubits (12 6 parameters ) 

5 qubits ( 62 parameters ) 
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FIG. 7. (Color online) Density distribution of the scaled fi- 
delity z for Bell and GHZ states. The measurement protocol 
is based on tetrahedron. Sample size n is equal to one million. 



with increasing the number of qubits. 



F=0,999 F=0.9999 




z=log(1-F) 

FIG. 6. (Color online) Density distribution of the scaled fi- 
delity z for three qubit state with different value of mixture 
GHZ and "white noise" states. The measurement protocol is 
based on dodecahedron. Sample size n is equal to one million. 

And finally, Fig. [7] presents the behavior of the recon- 
struction precision for Bell and GHZ states when the 
number of qubits grows from two to eight and the sample 
size is one million. The measurement protocol is based 
on tetrahedron. It is evident that the precision of recon- 
struction falls as rapidly as the width of the distribution 



V. EXPERIMENT AND ANALYSIS 

For testing the theoretic approach, described in the 
previous sections, we have prepared a whole family of 
polarization states for qubits and ququarts. In particular 
we considered both pure and mixed qubit states and pure, 
mixed, entangled and separable states of ququarts. 

Qubits. The experimental setup for generation and 
measurement of qubit states is shown in Fig. [5] Set- 
ting up the experiment we pursuit the following goals: 
(1) The setup should allow performing the transition be- 
tween pure and mixed polarization states of qubits. (2) 
Three different measurement protocols (R4, K4, B36) 
should be realized. (3) The setup should allow perform- 
ing measurements with different sample size (statistics). 

To satisfy the first request we used a wide-band 
light source (incandescence lamp) passing through a 
monochromator and thick birefringent plates. This pro- 
vides a variable spectral range around the chosen central 
wavelength 1.55 /im and a controllable phase delay be- 
tween basic polarization modes. Therefore, by changing 
the spectral range of the light (within 1-23 nm) we were 
able to vary the purity of the output polarization state. A 
parallel light beam was formed by a SMF28 single-mode 
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fiber with F240FC-1550 micro-objectives placed at its 
input and output. A Glan-Thompson prism was used for 
preparation of the horizontally polarized state \H) serv- 
ing as initial state for the following transformations. As 
a result, the original pure state was transformed into a 
mixed state with a degree of purity, depending on the 
spectral width of the detected radiation. A smooth tran- 
sition from pure states to total mixtures was achieved 
with increasing spectral width of radiation. Since the 
"measurement part" of the setup had a finite spectral 
band, the polarization states at different spectral com- 
ponents within this band were integrated, which corre- 
sponded to the registration of mixed-polarization state. 



H 45" 



lamp 



monochromator 



Preparation part 




FIG. 8. (Color online) Experimental setup for different to- 
mographic reconstructions of qubits with variable degree of 
mixture. Here H, V are Glan-Thompson prisms, ho is a thick 
quartz plate, hi is a thick tilt plate, hd is a driving phase 
plate, /i2,3 are tomographic phase plates, and D is a detector. 

For the preparation of pure polarization states of 
qubits we selected a spectral range about 1 nm. The ini- 
tial state \H) passed through a phase plate hd (436 /im) 
oriented at 0°, 45°. So the states under consideration had 
the following form respectively: 



= \H), |* 2 ) = 0.345|if) - 0.939i|F). 



(44) 



The second requirement is achieved with a standard 
quantum tomography method that used two achromatic 
quartz plates /i2=441 /im, /i3=313 //m served for recon- 
struction of polarization states. These plates were ori- 
ented at particular angles a, (3 (with respect to the verti- 
cal axis), so light passing through the plates and follow- 
ing vertical polarizer was projected onto the necessary 
set of states required by protocols J4 and R4. In the first 
protocol J4, suggested in [f|, projective measurements 
upon some components of Stokes vector were performed: 
\H), \V), |45°) = j=(\H) + \V)), \R) = ±(\H) + i\V)). 
Basically the projective measurements can be chosen ar- 
bitrarily, in particular if the measured qubits were pro- 
jected on the states possessing tetrahedral symmetry 
then the protocol transforms to R4. There are several 
works showing that due to the high symmetry such pro- 
tocol provides a better quality of reconstruction |23j-|25| . 
Another protocol B9 exploits a single plate and fixed po- 
larizer. The corresponding measurements have been per- 
formed for each of the nine orientations of the plate with 
a step of 20°. We have chosen the "optimal" thickness of 
the plate /i3=313 /jm and achieved the better condition 
number (K ~ 2.7) for reconstruction of mixed states. 
Also for testing theoretical predictions with this protocol 
we used other plates /i3=824 /im and /i3=358 /zm and 



corresponding condition numbers were -Ks24 = 14.9 and 
K 358 = 183.9. 

The third requirement was satisfied by using a single- 
photon detector based on GaAs-based avalanche photo- 
diode with a pigtailed input and an internal gate shapcr 
[I?} . The total number of registered events was varied by 
changing the gate rate at a fixed gate width (about 10 ns) . 
This simply follows from the fact that the fixed photon 
flux comes continuously and increasing (decreasing) the 
gate width just increases (decreases) the probability to 
register a photon. 

The goal of our first experiment was the optimiza- 
tion of a quantum tomography protocol for polarization 
qubits. As an example we selected protocol B9 imple- 
mented by means of a single phase plate and a fixed 
polarizer. Note that the projection in the protocol B9 
is performed onto non-fixed states, in contrast to the 
protocols J4, K4, and depends on the parameters of the 
plate, that is, on its optical thickness, 5 = irhAn/X, and 
the orientation angle /3, where An is the birefringence 
of the plate material at a given wavelength A and h is 
its geometric thickness. It is worth to mention that the 
parameters of this protocol can be chosen by doing the 
experiment, depending on available resources: in some 
sense this choice can be done in an optimal way, but us- 
ing the plate with optimal condition number. Figure [5] 
shows the calculated condition numbers K and maxi- 
mum losses L on the Poincare-Bloch sphere as functions 
of the optical thickness of the phase plate for the B9 
protocol. The dependencies arc periodic with a period 
of 7r. Poor conditionality occurs at 8 — 0,-|,7r. Both 
quantities K and L tend to infinity at these particu- 
lar points. For this protocol the best (lowest) achiev- 
able parameter K takes the value 1.85. For example, to 
reach such a value one might choose plates with optical 
thickness 5 = 0.356-7T or 8 = 0.6447T. The corresponding 
condition numbers for protocols J4, R4 are Kja 3.23 
A'r,4 = \/?> ~ 1.73. Particularly for our experiment we 
have chosen a set of following three plates: h^—313 /im 
(optimal), /i3=824 /im (medium), h^=358 /im (nonop- 
timal) with condition numbers K = 2.7, K = 14.9 
and K = 183.9. The optimal value of maximum losses 
L op t = 1-47 occurs at the points S = 0.3917T, 0.6097T. The 
corresponding losses are somewhat lower than those for 
the R4 protocol (L = 1.5), and essentially lower than 
those for the J4 protocol (L = 4.52). Note that these 
curves on Fig. [9] are similar in form, but their minima 
correspond to slightly different values of 6. However, the 
singular values of K and L correspond to the same 5 
value. Thus, both criteria allow separating the regions of 
the protocol parameters for which the results would be 
unsatisfactory. Unfortunately, for a large dimension of 
the Hilbert space for reconstructed states, the scanning 
of L values is a complicated computational problem. For 
this reason, to optimize such protocols, we suggest to use 
the parameter K . 

Figure [TU1 demonstrates the crucial difference between 
the optimal and nonoptimal approach for the exam- 
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FIG. 9. (Color online) Calculated condition numbers K and 
maximum losses max(L) on the Poincare-Bloch sphere versus 
the optical thickness 8 of a phase plate for the B9 protocol. 
The dashed and dash-dotted lines are for the R4 and J4 pro- 
tocols, respectively. 



pie of considered protocol B9. Here the optimal pro- 
tocol [Figs. ITDT a) and HUT b)] corresponds to the plate 
with thickness /i=293.8 /im and condition number K = 
2.2548. This protocol provides the minimum of maxi- 
mally possible losses (minimax). The nonoptimal proto- 
col [Figs. [TUtc) andlTOTd)] corresponds to the plate with 
thickness /i=358 /im and condition number K — 183.9. 
The sample size is 10 4 . Figures ITOT a) and fTOlf b) . which 
correspond to the optimal protocol, describe the case of 
guaranteed number of nines in Fidelity being not less 
than z m j n = 3.83 (so that the mean fidelity (f22j) for 
all states on Bloch sphere can not be lower). At the 
same time, Figs. ITOT c) and ITOT d). which correspond to 
the nonoptimal protocol, stand for a very low level of 
fidelity (z min = 0.065), which means that almost none 
of the states on Bloch sphere can be reconstructed us- 
ing the protocol. Figures fTDTb) and [TUTd) are different 
from Figs. [TUT a) and [TUtc) by the step of the plate ori- 
entation [20 degrees for Fig fToT a) and lTOT c) and 1 degree 
for Figs. IT07 b) andfTOld). respectively], so that the num- 
ber of rows in the protocol grows from 9 to 180. We see 
that Fig. [TOtb) defines a more symmetric distribution of 
precision compared to Fig. [TUta). even though the range 
of values [z m i n , z max ] barely changes. A comparison of 
Fig. [TUtc) with Fig. ITOT d) shows that in the nonoptimal 
case increasing the number of projections little affects the 
distribution of Fidelity on Bloch sphere. 

For protocol B9 the statistical reconstruction of the 
prepared states f|44[) has been performed at given sam- 
ple sizes. One of the considered states belongs to 
an area of small losses for all three plates, and the sec- 
ond state 1^2) was changed out of this area. As an 
example, we present Fig. [TTJ that shows the calculated 
widths of fidelity distributions at 1% and 99% quantilcs 
as well as the experimentally reconstructed values for the 
states (|4"4")) . Two polarization states ([4"4"f were measured 
and reconstructed using three phase plates chosen above. 
Experiments were performed with different sample sizes; 




0.5 1 1.5 2 0.5 1 1.5 2 




FIG. 10. (Color online) Difference between the (a) and (b) 
optimal and (c) and (d) nonoptimal approach for protocol B9. 
(a) and (c) and (b) and (d) correspond to 20° and 1° steps of 
the plate orientation, respectively. 



that is, the total number of the pulses coming from single- 
photon detector in a fixed time was varied and served as 
a parameter of the problem. 

The approach described above provides the ideal ac- 
curacy level for quantum state reconstruction. It means 
that the fluctuations of the estimated quantum states 
cannot lead to uncertainties smaller than this limit. The 
presence of instrumental errors and uncertainties makes 
this level to be exceeded. Indeed, Fig. [IT] shows that, 
above some sample size, the experimental value of fidelity 
falls out the theoretical uncertainty boundary shown 
as dotted lines (corresponding to 1% significance level 
which characterizes a given protocol). This happens 
since instrumental uncertainties prevail over the statis- 
tical ones and indicates that either state preparation 
stage or measurement procedure were not performed ac- 
curately enough. It is clearly seen that both theoretical 
distributions and experimental points for protocol with 
nonoptimal choice of plates strongly depend on an initial 
state. The quality of reconstruction (fidelity) is varied in 
the range 0.6030-0.9999 for both reconstructed states for 
the protocol with the plate /i3=358 ^m (nonoptimal) and 
in the range 0.9950-0.9997 for the plate /i 3 =313 fxm (op- 
timal) . It means that if the set of initial states is known 
(as it happens in quantum process tomography) one can 
choose the parameters of the protocol that are not opti- 
mal for an arbitrary state, but provide the highest accu- 
racy for the selected set of states. In both cases, when 
the input states are unknown, like in quantum state to- 
mography, it is better to use optimal choice of the plates 
with minimal condition number, that will ensure high 
accuracy of the reconstruction of arbitrary states. 

Mixed States: modeling and experiment. The goal of 
the second experiment with qubits was the generalization 
of the developed approach to the family of mixed states of 
polarization qubits [13 ■ As an example we have tested 
three protocols again, namely J4, R4, and B36. The 
protocol B36 has been considered at optimal parameters 
(phase plate /i3=313 /mi). This protocol is similar to 
the one described above (B9), but measurements were 
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FIG. 11. (Color online) Reconstruction of pure qubit states by protocol B9 at different thicknesses of the plate. Vertical bars 
show 1% and 99% quantiles for fidelity distributions. Dotted lines connecting lower bar ends point out critical significant levels. 



performed at 36 consistent orientations of the phase plate 
(0°-360°, step 10°) instead of 9. 

As we have pointed out before, the condition num- 
bers K for R4 and J4 protocols take the following values: 
Kr4 = V3 « 1.73, K.n m 3.23 while for B36 it becomes 
-Kb36 ~ 2.7. Thus, we expect that the symmetrical proto- 
col i?4 provides with better state reconstruction quality. 
We have checked this statement with numerical simula- 
tions of each protocol applied to qubit states with various 
degree of mixture depending on the sample size. 

For the preparation of mixed states we started from a 
pure state = \H) which passed through one or two 
thick bircfringent plates, oriented at 45°, where partial 
decoherence between vertical and horizontal basic polar- 
izations took place. For that purpose we used a quartz 
plate with a thickness of ho=10 mm and additional tiff 
plate with a thickness of h\=A mm. Changing the width 
of the spectrum it is possible to pass from completely 
coherent (pure) case to a very narrow spectral line to a 
completely incoherent (i.e. mixed in polarization), when 
the difference between the optical lengths for two orthog- 
onal polarizations exceeds the coherence length of the 
radiation under consideration: 

h (n± - n||) quartz + hi(n± - n\\) tis ) > l coh , (45) 

where 

A 2 

hah ~ 2^ ~ 100 A ™' 

Thus, at the output of plates ho, hi components with 
orthogonal polarizations evolve with a random phase, 
which leads to a mixed state. 

For determining the accuracy of quantum tomography 
and fidelity, we need to calculate the density matrix of 
a prepared mixed state. For this purpose we divide the 
frequency spectrum of the radiation in small parts and 



represent the polarization state of a qubit as a superpo- 
sition of states corresponding to different frequencies in 
the spectrum: 

l*> = EjkOfc|*(wfc)>, m 

\<f>(cj k )) = ci(uj k )\H)+c 2 (u; k )\V}, 

where amplitudes a k are defined by the spectrum shape. 
The density matrix of the state before transformation 
(gBJ has the form: 

r = |*)<*| =X)o fc o5|$(o> fc )X*(w j )| (47) 

The next step of modeling includes the calculation of 
phase plates action on the qubit state. The unitary trans- 
formation on state (|47p is given by matrix 

Gk) = f *\ " fe ] (48) 
V - r k h J 

where 

t k = cos 5 k + i sin 5 k cos 2a, 
r k = isin5 fc sin2a, S k = ^(n^ - n^)h/X k . 

Here t k and r k are the amplitude transmission and re- 
flection coefficients of the wave plate at fixed frequency, 
5 k is its optical thickness, h is the geometrical thickness, 
a is the orientation angle between the optical axis of 
the phase plate and vertical direction. The measurement 
part of the experimental setup does not distinguish fre- 
quency modes, so the theoretical polarization mixed state 
will be described by a reduced density matrix: 

p = ]T \a k \ 2 G(u k )\$(u k )){$(uj k )\G + (Lj k ). (50) 

k 
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FIG. 12. (Color online) Theoretical distribution of the real and imaginary parts of the density matrix for states with varying 
degrees of purity. The spectral widths of the spectrum are (a) 1.6 nm, (b) 7 nm, (c) 13 nm. Distribution (d), for 22 nm, 
corresponds almost to a completely mixed state. 



The corresponding density matrixes are obtained by 
integrating formula (|50|) with the distribution of weights 



\ak\ as a function sine (x). Figure 1121 shows the graph- 
ical representations of the real and imaginary parts of 
the theoretical qubit density matrices at different widths 
of the spectrum: 1.6 nm, 7 nm, 13 nm, 22 nm. When 
increasing the width of the spectrum the non-diagonal 
components, responsible for the correlation, decay and 
the state becomes completely mixed. 

In the experiment we have prepared four polarization 
states of qubit with following degrees of mixture: 3% 
(1.6 nm), 30% (7 nm), 66% (13 nm), 100% (22 nm), 
where the state purity is analyzed by calculating the state 
entropy defined as S = — 5Z n =l ^ n 1°§2 ^n, A„ being the 
eigenvalues of the density matrix p. For each protocol 
various measurements at different sample sizes were per- 
formed. The results are shown in Fig. Q2J 

Figure [T3] presents the theoretical distributions of fi- 
delity at 1% and 99% quantiles for each considered pro- 
tocol. Experimental values are indicated by points. The 
dashed lines connecting the 1% quantiles show a theoreti- 
cal lower boundary of the fidelity distribution and depend 
on the type of protocol: the lower line, the greater con- 
dition number and the function of losses. Fig. [13] demon- 
strates that for a small sample size statistical uncertain- 
ties dominate over the instrumental ones. Starting from 
sample size about (2—3) x 10 3 the experimental points fall 
out the theoretical uncertainty boundary and the proto- 
col achieves a coherent sample size. This means that in- 
strumental uncertainties (the accuracy of setting angles, 
thicknesses of phase plates, etc.) prevail over the level 
of statistical errors. From this point a further increase 
of sample size does not improve the quality of the recon- 
struction of the quantum state. However, a comparison 
of experimental results with the theoretical distribution 
serves as effective method for setup adjusting, stability 
of the preparation/measurement systems, etc. 



Furthermore, from Fig. 1131 one can evince that the ac- 
curacy of the reconstruction increases with the degree of 
mixture: for the states close to a pure ones (3% mixture) 
the fidelity achieves an average value 0.93-0.94, while 
the best accuracy is achieved for totaly mixed state cor- 
responding to the center of the Bloch sphere, F > 0.99. 
The theoretical estimates are consistent with both nu- 
merical and real experiments. 

Ququarts: pure and mixed states, preparation and mea- 
surement. For a further verification of our approach we 
prepared also a family of quqarts biphoton polarization 
states [i^ , which can be easily converted into either en- 
tangled (in polarization) or factorized states, both pure 
and mixed. 



p = (l-p)\% v 



pure | 



H X H 2 YH 1 H 2 \ 



where 



\V X V 2 ){V 1 V 2 \), 



c 1 \H 1 H 2 ) + c 2 e^\V 1 V 2 [ 



(51) 



(52) 



with real amplitudes c\ and c 2 and relative phase shift 
if. The coefficient p defines the degree of mixture. 

The set-up is schematically depicted in Fig. [TU For 
generation of pure biphoton-based polarization ququarts 
(p = 0) we used a set of two orthogonally oriented 
type-I BBO crystals (1 mm), cut for collinear, frequency 
non-dcgcncratc phase-matching around the central wave- 
length of 702 nm. The crystals were pumped by a 
600 mW cw- argon laser operating at 351 nm. The Glan- 
Thompson prism (V) at vertical polarization and the 
half-wave plate (A p /2) placed in front of crystals allowed 
rotating the polarization of the pump by the angle 0, 
which controlled the real amplitudes c\ and c 2 in (|52p . 
A set of quartz plates QP introduced the relative phase 
shift tp between the horizontal and the vertical compo- 
nents of the UV pump. If <fi = we prepared the state 
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FIG. 13. (Color online)Reconstruction of qubits states with various degree of mixture. Vertical bars show 1% and 99% quantiles 
for fidelity distributions. Dotted lines connecting lower bar ends point out critical significant levels, (a) mixed 3%, (b) mixed 
30%, (c) mixed 66%, and (d) mixed 100%. 



I*) = \ViV 2 ), if = 22.5° and <p = 3tt/2 then the state 
transformed to = (\H X H 2 ) - |ViF 2 )) = I* - )- To 
maintain stable phase-matching conditions, BBO crys- 
tals and QP were placed in a closed box heated at fixed 
temperature. The lens L coupled SPDC light into the 
monochromator M (with 4 nm resolution), set to trans- 
mit "idler" photons at 710 nm. The conjugate "signal" 
wavelength 694 nm was selected automatically by means 
of the coincidence scheme. 

In order to prepare a mixed state it is necessary to 
introduce quantum distinguishability among the bipho- 
ton ququart basis states in a controllable manner. For 
the preparation of states with various degree of mix- 
ture the double-crystal scheme was complemented by a 
quartz plate with thickness /io, which was putted in the 
reflected arm of Brown- Twiss scheme. A thick quartz 
plate with vertically oriented optical axis introduced a de- 



lay between vertically and horizontally polarized photons 
that led to their temporal distinguishability, and hence, 
to the emergence of a mixture. Changing the width of 
the frequency spectrum, which, in our case, was deter- 
mined by the width of the monochromator slit, or the 
plate thickness it was possible to change the visibility of 
polarization interference. Our goal was the preparation 
and reconstruction (by several protocol, including non- 
optimal ones) of a set of states with different degree of 
entanglement and set of mixed two-qubit states. 

The reconstruction part for the ququart is shown in 
Fig. [TU First, the photon pair that forms the biphoton- 
ququart was split into two spatial modes by using a 
beam splitter BS. A photon in each arm, then, under- 
went a polarization state transformation with the use of 
a couple of zero-order wave plates (A/2, A/4) for proto- 
cols. Two Si-APD's linked to a coincidence scheme with 
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FIG. 14. (Color online) Experimental setup for different to- 
mographic reconstructions of photon pairs with variable po- 
larization entanglement and degree of mixture. Ar laser: ar- 
gon laser with wavelength 351 nm, M: mirror, V: vertical 
oriented Glan-Thompson prism, BBO: nonlinear Barium Bo- 
rate crystals, UVF: ultraviolet filter, A p /2, A p /4: half-wave, 
quarter-wave plates, L: lens with focus 20 cm, BS: beamsplit- 
ter, IF: interfilter, D: avalanche photodetectors, CC: coinci- 
dence circuit. 
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FIG. 15. (Color online) Dependence of \og(K) on plates thick- 
nesses for protocol B144. 



1.5 ns time window were used as single photon detectors. 
Registering the coincidence rate for different projections, 
that were realized by the half- and quarter-plates and a 
fixed analyzer, was possible to reconstruct the polariza- 
tion state of ququarts by protocols R16, J16, described 
above. This scheme also allowed reconstructing any ar- 
bitrary polarization ququart state by protocol B144. For 
the reconstruction of ququarts pure states by protocol 
B144 we used two retardation plates Wpl, Wp2 placed 
before beam splitter and deleted half-wave and quarter- 
wave plates. Unfortunately, this method can not be ap- 
plied for the reconstruction of mixed states of ququarts. 
This is due to the feature of preparation for mixed states: 
all projection measurements must be done after state 
preparation that is violated for such configuration of the 
tomography protocol. To solve this problem, we have 
used an equivalent scheme. In each arm we have es- 
tablished a couple of identical plates Wpl, Wp2 instead 
of half-wave, quarter-wave plates and rotated them syn- 
chronously. This scheme is completely equivalent to the 
standard protocol B144, but slightly more complicated 
in realization. By analogy to protocols B9 and B36 with 
single retardant plate we can choose the optimal param- 
eters for protocol B144. 

Let us analyze how the choice of plates thicknesses 
affects the quality of state reconstruction. 

Figure IT51 presents the calculated condition number K 
in common logarithmic scale for the set of thicknesses of 
the first and second quartz retardant plates. Specifically 
the first plate Wpl was varied within the limits from 
1.100 to 1.350 mm with a step 0.001 mm and a second 
plate, Wp2, was varied within the limits from 0.300 to 
0.550 mm with the same step. Here the thicknesses of the 
first and second quarts retardant plates are shown along 
the axis x, y, correspondingly. Blue color indicates areas 
with low value of condition number K, while the red one 
indicates areas with high value of K . It is obvious, that 



the arbitrary choice of plates most likely will be wrong 
from the point of view of completeness of matrix B and 
protocol optimality. In other words this figure serves as 
some sort of a "navigation map" for the measurement 
protocol. So to achieve good quality of the measurement 
one should select the plates thicknesses in blue areas and 
avoid red areas. For the present protocol the minimum 
achievable K is 11.35. For example, to reach such a value 
one can choose plates with thickness Wpl=1.207 mm and 
Wp2=0.520 mm. In our specific experimental configu- 
ration we have selected the plates Wpl=1.303 mm and 
Wp2=0.440 mm. Such thicknesses of the plates for the 
protocol B144 were chosen to show the difference in the 
quality of reconstruction. 

For all three protocols used in the experiment we have 
calculated the condition numbers K, which assume the 
following values: ivRi6 = 3, Kjiq ps 10, Kbim ~ 60. 
Thus, we expect that the symmetrical protocol R16 pro- 
vides with better state reconstruction quality (26j and the 
protocol B144 provides the worst one. We have checked 
this statement with numerical simulations of each pro- 
tocol applied to different two-qubit states depending on 
the sample size. 

As an example, let us consider the numerical recon- 
struction of the Bell state |$~) = ^= (|iTi-ff 2 ) - \ViV 2 )). 
Figure [If)] shows average fidelities, calculated according 
to (|2T|) as functions of a sample size for each protocol. It 
turns out that the difference between curves disappears 
at sufficiently large sample size (10 5 )[5(|, but for the same 
quality of state reconstruction the correct choice of pro- 
tocol allows using a smaller set of statistical data, i.e. fi- 
nally reduces total acquisition time. FigurefTolshows that 
protocols are ranged in accuracy as following: R16, J16, 
and B 144 in complete agreement with the range given by 
condition number K. Figure [T71 presents accuracy distri- 
butions calculated according to ([2~Tj) for the sample size 
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FIG. 16. (Color online) Dependence of the average fidelity on 
number of registered events forming the sample for Bell state 
reconstructed by protocols J16, K16 and B144. 



3 x 10 3 . 

It is clearly seen that the density distribution for the 
R16 protocol is narrower than the one for J16 and B144 
and localizes in the region of lower losses or higher fideli- 
ties. The distribution for B144 is broader and lower in 
comparison with R16 and J16. Obviously, the narrower 
distribution of fidelity indicates a better reconstruction 
quality in the sense that the reconstruction procedure 
returns a better-defined state. Thus, Figure [TTl confirms 
our expectation based on estimation of condition numbcr 
K; R16 achieves the best results. 

Let us consider now the experimental results. As an 
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FIG. 17. (Color online) Density distribution of the scaled 
fidelity z (lower abscissa) at n = 3 x 10 3 for Bell state re- 
constructed by protocols J16, K16 and B144. Upper abscissa 
presents regular fidelity. 



example, Fig. [18] shows both experimental points and 
theoretical distributions for the set of pure ququart states 



with various degree of entanglement: C, 



o, a 



¥>2 



0.83, C„ 4 = 1 correspondingly. Here C lfi is 



0.66, C V3 

quantity concurrence, which is defined as C = 2\c\c± — 
C2C3I and C = for a separable state and C — 1 for a 
maximally entangled state. 

Figure flUl presents experimental points and theoretical 
distribution for the set of mixed ququart states with vari- 
ous degree of mixture: p = 1, 0.3, 0.7 (p = 1 corresponds 
to totaly mixed state, p = 0.3 means that state contains 
30% of mixture.) 

On both figures the agreement between the experimen- 
tal results and a priori calculated expectations is note- 
worthy. 

Figure [T9l shows that the accuracy of reconstruction in- 
creases with the degree of mixture. However from a sta- 
tistical point of view there is a problem. A small weight of 
the mixed states, adding to the original pure state, leads 
to a small change in the statistical data while the number 
of parameters for ququart state increases from 6 to 15. 
As a result the amount of new information that arises 
in quantum measurements is not sufficient to achieve an 
adequate estimation of all the parameters. 



VI. DISCUSSION 

Many remarkable works in quantum state tomogra- 
phy developed different theoretical approaches. In the 
present section we shall briefly discuss the features and 
advantages of our approach in this framework. 

One of the main issues in quantum tomography is the 
choice of adequate parameterization for quantum states. 
Bloch sphere representation is one option. In this repre- 
sentation density matrix p is defined in Hilbert space of 
dimension s by the following equation [H, l35l . l5lT | 



P 



l -E 

s 



(53) 



where E is the s x s identity matrix, Uj is a vector in 
real Euclid space of size s 2 — 1 , and j changes from 1 to 
s — 1. Respective Hermitian basis matrices <7j meet the 
following conditions: 



tr(crj) = 0, ti(<jjO- k ) = S jk . 



(54) 



Note, that for singlc-qubit states when s — 2, matrices Oj 
coincide with ordinary Pauli matrices (multiplied by 77^)- 
The radius of 3D Bloch sphere in this representation is 
also equal to 

One may easily calculate Bloch vectors Vj from the 
known density matrix p: 



tT(po-j). 



(55) 



However, Bloch parameterization has a significant draw- 
back. For every real vector v Eq. (|55|) guarantees that 
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FIG. 18. (Color online) Reconstruction of ququart states with various degree of entanglement. Vertical bars show 1% and 99% 
quantiles for fidelity distributions. Dotted lines connecting lower bar ends point out critical significant levels. 



the density matrix p is Hermitian, but does not guarantee 
that it is positively defined. As a result, we can obtain a 
non-physical density matrix using Eq. (|53[) . Separation 
of admissible density matrices from those that can not 
physically exist is a non-trivial task for all dimensions 
higher than two. In particular, for s > 2 Bloch repre- 
sentation is practically useless for description of matrices 
of non-full rank, for which the number of degrees of free- 
dom of quantum state (|36j) is less than dimension of Bloch 
vector s 2 — 1. This drawback of Bloch parameterization 
complicates the use of Fisher information matrix and cal- 
culation of statistical estimates precision in general. 

An adequate and convenient procedure of parameteri- 
zation is important for numerical procedures of calculat- 
ing density matrix from experimental data because for 
every iteration the approximate solution must lie in the 
range of physical density matrices only. For instance, 
such issues may arise for problems of estimation of quan- 
tum state by maximum likelihood method. 

The problem of estimating a density matrix by max- 



imum likelihood method was considered in 

The respective likelihood equation has the following form 

Rp = p, (56) 

where R is some operator. 

The iteration procedure from jth to j + lth step for this 
equation has the following form: p> +1 = Rp 3 . Such pro- 
cedure, however, does not even provide Hermitian prop- 
erty for density matrix. 

Another procedure that is widely used [5l[ pi +l = 
\iR-P 1 + p' J ' R) provides the Hermitian property but does 
not guarantee positive definiteness in general. Fortu- 
nately there are known robust, well-converging proce- 
dures providing with both Hermiticity and positive defi- 
niteness simultaneously {EH ]. 

In this paper we consider a natural approach for quan- 
tum theory that is based on re pre senting the quantum 
state as a mix of pure states [l7l Il8l. |40| . Consider trans- 
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FIG. 19. (Color online) Reconstruction of ququart states with various degree of mixture. Vertical bars show 1% and 99% 
quantiles for fidelity distributions. Dotted lines connecting lower bar ends point out critical significant levels, (a) p = 0.3, (b) 
p — 0.7, and (c) p = 1. 



formation of the density matrix to diagonal form. 

p = UDU^. (57) 

Here U is an unitary matrix (its columns define eigen- 
vectors of density matrix) , D- is a diagonal non- negative 
matrix (its diagonal is formed by eigenvalues of density 
matrix that we shall put in decreasing order). Let the 
rank of density matrix be r (1 < r < s). We shall elim- 
inate all non-zero rows and columns of matrix D and 
shall leave only the first r columns in matrix U. Then 
the density matrix could be represented in the following 
compact form. 

p=LL\ (58) 

where L = UDi . 

Here L is a complex matrix of dimension s x r that 
defines the actual parameterization of density matrix. 
It represents a purified state vector (probability ampli- 
tudes). For representation as a column vector, one should 
simply place the second column under the first one and so 
on. Note that the purified state is defined ambiguously, 
because the density matrix does not change during trans- 
formation 

£->•£' = LV, (59) 

where V- is an arbitrary unitary matrix of size r x r. 

It is crucial that the ambiguity (|5"9"|) does not affect the 
procedure of statistical reconstruction of density matrix, 
because all possible purified state vectors define the same 
density matrix. 

We can choose the unitary transformation V such that 
in matrix L' all elements on the principal diagonal will 
have real strictly positive values, while elements higher 
the principal diagonal shall be equal to zero. 

For a density matrix of full rank, such matrix L' shall 
form a lower triangular matrix. In mathematics such 
decomposition is usually called Cholcsky decomposition. 
Representation of density matrix in such form was origi- 
nally studied in [l(| ■ Note that from the point of view of 



physics, Cholcsky representation is just one of the possi- 
ble ways of recording a purified quantum state. 

A purified state vector represents the most simple and 
most natural in view of physics form of parameterization 
of density matrix. It is important that parameterization 
based on purification of quantum state radically simpli- 
fies the theory of statistical estimates. 

For example, the precision of estimates by maximum 
likelihood method is defined by matrix of full information 
that is analogous to Fishers information matrix in rela- 
tion to estimation of quantum state vector [13, H3, H(j ■ 

The matrix of full information is defined by the follow- 
ing equation: 

g = 2E MMV)!, (60 ) 

3 J 

The object is defined in real Euclidian space of dimension 
2rs. To obtain state vector c in this representation, one 
should place the imaginary part of purified state vector 
under its real part. Intensity matrices Aj represent ex- 
pansion of matrices © to the considered Euclid space. 
The sum in Eq. (|60[) is taken for all measurement proto- 
col rows. Matrix H becomes a real symmetric matrix of 
dimension 2rs x 2rs. 

If the quantum measurement protocol is complete, 
then for the information matrix H vjj = (2s — r)r out 
of Its eigenvalues are strictly positive, while the other 
r 2 are equal to zero. One may formulate a universal sta- 
tistical distribution (|21j) based on information matrix H 
that will describe precision of statistical reconstruction 
of quantum states. 

Eigenvectors of matrix H corresponding to non-zero 
eigenvalues define directions of fluctuations of a quantum 
state and its norm, while eigenvectors that correspond to 
zero eigenvalues define directions of insignificant fluctua- 
tions that are due to the ambiguity ([51?]) in definition of 
the purified state vector. 

While information matrix H defines detailed charac- 
teristics of precision of reconstruction of an arbitrary 
quantum state, the measurement matrix B defines 
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the quality of the protocol as a whole. The condition 
number K separates complete and well defined quantum 
measurement protocols from incomplete and ill-defined 
ones. 

The approach developed above is based on analysis of 
completeness, adequacy and precision of quantum mea- 
surements and it could be applied to analysis of quality 
and efficiency of arbitrary protocols of quantum measure- 
ments. 

As a simple example, let us consider a single-qubit pro- 
tocol pro pos ed in works [HI, HU (see Fig. 1 in [35| and 
Fig. 3 in [53j). The considered approach is an optimized 
measurement of state 



We see that for the state 




Measurements are conducted by means of a set of half- 
and quarter- wave plates that lead to projections to the 
following basis states 




the protocol is indeed the optimal one because the losses 
are at minimum L m j n = 1. However for a number of other 
states the protocol is rather not optimal (maximum losses 
are quite large L max ~ 16.84). For instance, the protocol 
considered here has much less precision than those stud- 
ied in section IVT1 which are based on symmetries of cube 
and octahedron and which also have few rows (6 and 8 
respectively), but provide much lower precision losses in 
narrow range from L m { n = 1 to L max = 9/8. 

Note that the condition number K for the protocol is 
approximately 4.07 times higher than for protocols based 
on regular polyhedrons, which also certifies its lower qual- 
ity. 



VII. CONCLUSIONS 
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where h and q are orientation angles for half-wave and 
quarter-wave plates respectively. This measurement de- 
fines two strings of instrumental matrix X . The authors 
of the protocol proposed a set of four measurements of 
such type (8 rows in total). Measurements are defined by 
the following orientation angles for the plates: hi =0°, 
qi = 0°; h 2 = 20°, q 2 = 45°; h 3 = 25°, q 3 = 45°; 
h 4 = 45°, q 4 = 0°. 

A complete analysis of precision losses (fM)) during re- 
construction of single-qubit pure quantum states using 
the protocol is given in Fig. [2D] (little precision loss states 
are marked blue, large precision losses a marked brown). 
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FIG. 20. (Color online) Analysis of precision losses for proto- 
col proposed in [H, [k| ■ 



In this paper, we proposed a methodology to assess the 
quality and efficiency of quantum measurement proto- 
cols. It was shown that the proposed approach, based on 
the analysis of completeness, adequacy and accuracy, can 
be successfully applied to arbitrary quantum states and 
measurement protocols. Analysis of the completeness of 
the protocol allows to answer the question: Is the set 
of operators considered sufficient to assess an arbitrary 
pure or mixed quantum state? The operational crite- 
rion of the protocols' quality follows from the analysis 
of completeness. The criterion is based on the condition 
number of the measurements matrix K. An analysis of 
the adequacy allows one to answer the question about the 
correspondence between the statistical experimental data 
and the quantum-mechanical mathematical model. The 
analysis of the adequacy guarantees the correctness of the 
following procedures of statistical quantum state recon- 
struction by using the maximum likelihood method. The 
developed method of statistical reconstruction is based 
on on the procedure of quantum state parameterization 
using the procedure of purification. Such parametriza- 
tion allows us to introduce a universal distribution for Fi- 
delity, which provides all information about the accuracy 
of the the quantum state reconstruction. Various numer- 
ical examples and results of physical experiments have 
been analyzed for testing theoretical predictions both for 
pure and mixed states and demonstrating approach va- 
lidity. 
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